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Background: Average energies of nuclear collective modes may be efficiently and accurately 
computed using a non-relativistic constrained approach without reliance on a random phase ap- 
proximation (RPA). Purpose: To extend the constrained approach to the relativistic domain and 
to establish its impact on the calibration of energy density functionals. Methods: Relativistic RPA 
calculations are compared against the predictions of the corresponding constrained approach using 
two accurately calibrated energy density functionals. Results: We find excellent agreement — at 
the 2% level or better — between the predictions of the relativistic RPA and the corresponding con- 
strained approach for magic (or semi-magic) nuclei ranging from ^^O to ^"^Pb. Conclusions: An 

C^ efficient and accurate method is proposed for incorporating nuclear collective excitations into the 

^~~^ calibration of future energy density functionals. 
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00 I. INTRODUCTION 



^S^ A tractable microscopic theory that both predicts and provides well-quantified theoretical uncertainties for nuclear 

I properties throughout the nuclear chart is the guiding principle in the quest of a Universal Nuclear Energy Density 

' fj Functional pQ [2] . Density functional theory (DFT) provides a powerful — and perhaps unique — framework for the 

^ accurate calculation of the ground-state properties and collective excitations of medium-to-heavy nuclei. Based on 

^ the seminal work by Kohn and collaborators [MS], DFT shifts the focus from the complicated many-body wave 

'~~' function to the much simpler one-body density. By doing so, the formidable challenge of deducing the exact ground- 

I state energy and one-body density from the many-body wave function "reduces" to the minimization of a suitable 

^ functional of the density. Of course, this enormous simplification comes at a price. Whereas DFT establishes the 

T-H existence of a well-defined density functional, it offers no guidance on how to build it. Nevertheless, by strictly 

CNl focusing on measurable quantities, DFT offers the opportunity for using physical intuition and symmetries/constraints 

\l to guide the construction of the functional. Moreover, Kohn and Sham have shown how the ground-state energy and 

^^ corresponding one-body density may be obtained from a variational problem that reduces to a self-consistent solution 

■^ of a set of mean-field-like {"Kohn- Sham") equations [J. However, although the form of the Kohn-Sham equations 

^^ is reminiscent of the self-consistent Hartree problem in the presence of an underlying (bare) nucleon-nucleon (NN) 

^^ interaction, the constants that parametrize the Kohn-Sham potential are directly fitted to many-body properties (such 

, . as masses and charge radii) rather than two-body data. In this manner the complicated many-body dynamics gets 

^ implicitly encoded in the parameters of the model. In principle, a proper implementation of DFT and the Kohn-Sham 

equations incorporates all many-body effects [5 . 

The construction of an energy density functional (EDF) starts with the selection of terms that incorporate important 
J^ symmetries and features of the nuclear dynamics. Although in principle such terms may be inspired by the underlying 
nucleon-nucleon dynamics, the explicit value of the coefficients in front of these terms (i.e., the model parameters) is 
customarily obtained through the minimization of a quality (e.g., x^) function. Thus, the model parameters may — and 
in general do — differ significantly from the "bare" NN values. This is a reflection of the fact that the parameters of 
the model are calibrated to physical observables that incorporate few- and many-body correlations. In the language 
of DFT, most of the complicated many-body dynamics gets implicitly encoded in the parameters of the model. Once 
the EDF has been properly defined, a set of accurately-measured ground-state observables is selected to constrain 
the model parameters through the minimization of a quality function. However, given that ground-state observables 
are fairly insensitive to fluctuations around the equilibrium density, one must often resort to "pseudo-data" — in the 



X 



'Electronic address: wc09c@my.fsu.edu| 
t Electronic address: .ipiekarewicz@fsu.edu| 



-f Electronic address: mariocent ellesQu b . edu 



form of various bulk properties of infinite nuclear matter — to better constrain the functional. Regardless, once the 
functional is properly defined and a fitting protocol established, minimization of the quality function defines the 
model. Until recently, once the minimum was found, one proceeded to validate the model against observables not 
included in the quality function [B]. Lately, however, a few studies have been devoted to map the landscape around 
the x^-minimum [7, .9,. Among the wealth of information revealed by such detailed statistical studies is the degree of 
correlation among various observables. Moreover, one can also use such a covariance analysis to assess the robustness 
of the model. For example, one could ask whether certain linear combinations of model parameters remain poorly 
constrained by the choice of observables. We find this to be particularly true in the case of the isovector sector that 
is hindered by the unavailability of highly accurate data on neutron skins [Ij . 

It is the main goal of the present contribution to explore the feasibility of supplementing ground-state observables 
with nuclear excitations in the calibration of the quality function of relativistic functionals. We are confident that 
such an approach will relax the need for pseudo-data while providing better constraints on the determination of the 
model parameters. Specifically, we advocate supplementing the quality function with centroid energies of monopole 
resonances for nuclei with a wide range of neutron-proton asymmetries. By doing so, the quality function becomes 
highly sensitive to the incompressibility of neutron-rich matter — which itself depends on both the compression modulus 
of symmetric matter and the slope of the symmetry energy |10H12j . 

Perhaps the most serious impediment to the implementation of this program is numerical in nature. By itself, 
computing the distribution of monopole strength for a wide range of excitation energies for models with finite-range 
interactions is numerically intensive [13]. Thus, embedding such RPA calculation into a complex x^-minimization 
routine becomes impractical even for today's most powerful computers. However, the computation of centroid energies 
relies on knowledge of just a few moments of the distribution. In particular, an estimate of the centroid energy may 
be obtained from knowledge of two moments: the energy weighted sum (rrii) and the inverse energy weighted sum 
(m_i). Remarkably, the fully self-consistent m_i moment may be calculated quite generally from the ground-state 
properties of a slightly modified (i.e., "constrained") density functional [I4| . Similarly, the mi moment may also be 
obtained by computing the ground-state expectation value of a suitable operator ((r^) in the case of the isoscalar 
monopole mode) 15J. Indeed, this constrained approach has been used successfully to test self-consistency violations 
in non-relativistic Skyrme-RPA calculations [16j . 

The constrained approach has been widely employed in the non-relativistic mean-field theory as a surrogate for 
explicit RPA calculations of giant monopole energies, see for example Refs. [T6ti23] and references quoted therein. In 
contrast, no more than a handful of constrained calculations have been reported in the context of relativistic mean-field 
theory. For example, relativistic constrained calculations of giant monopole energies of finite nuclei were performed 
in Refs. PH [2S] within the linear Walecka model. Some time later, constrained monopole energies using non-linear 
RMF models were obtained in Ref. [26], whereas a somewhat different constrained approach, the so-called generator 
coordinate method (GCM), has been applied in RMF studies of giant monopole resonances in Refs. pTl - BU] , The 
available relativistic literature with quantitative analyses of the degree of agreement between constrained calculations 
and the corresponding RPA results is scarce — we are aware of only the RMF constrained- vs- RPA study of the m^i 
moment of the monopole mode carried out in Ref. [W for ^"^Pb. On the other hand, to our knowledge, the constrained 
approach has not been systematically implemented in the context of relativistic mean-field theories with the goal of 
supplementing the calibration of nuclear functionals with collective excitations. The reason for such an imbalance in 
comparison with the non-relativistic framework may be fairly evident. Whereas the derivation of the inverse energy 
weighted sum rule hinges exclusively on the validity of perturbation theory, the derivation of the energy weighted 
sum rule (EWSR) assumes a non-relativistic kinetic energy operator (i.e., quadratic in the momentum). However, 
extending the constrained approach to the relativistic domain is plausible against the backdrop of density functional 
theory. DFT prescribes how the linear response of the system may be extracted from computing small density 
fluctuations around the variational minimum [5]. In particular, DFT makes no demands on whether the functional 
should be of a relativistic or non-relativistic character; the precise form of the functional then becomes a matter of 
convenience. 

The manuscript has been organized as follows. In Sec. |ll] we review the necessary formalism required to compute 
the excitation energy from the complete distribution of isoscalar monopole strength. In addition, we describe how 



the r7i_i and nii moments are computed in the constrained approach. In Sec. Ill we compare results obtained from 
the constrained approach against those extracted from the full distribution of strength. We conclude in Sec. IV with 
a summary of our results and plans for the future. 



II. FORMALISM 

The starting point for the calculation of the nuclear response is the interacting Lagrangian density of Ref. [31] 
supplemented by an isoscalar-isovector term first introduced in Ref. |32]. That is, 
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The Lagrangian density includes an isodoublet nucleon field {ip) interacting via the exchange of two isoscalar mesons, 
a scalar ((/>) and a vector (y^), one isovector meson (6^), and the photon {A^^) 031 IMl- In addition to meson-nucleon 
interactions, the Lagrangian density is supplemented by four nonlinear meson interactions with coupling constants 
denoted by k. A, C,, and Av. The first two terms {k, and A) are responsible for a softening of the equation of state of 
symmetric nuclear matter at normal density that results in a significant reduction of the incompressibility coefficient 
relative to that of the original Walecka model [331 133 [3S]- Indeed, such a softening is demanded by the measured 
distribution of isoscalar monopole strength in medium to heavy nuclei |131 1171 [371 I38j . Further, omega- meson self- 
interactions, as described by the parameter C, serve to soften the equation of state of symmetric nuclear matter at high 
densities and at present can only be meaningfully constrained by the limiting masses of neutron stars [35] . Finally, 
the parameter Av induces isoscalar-isovector mixing and is responsible for modifying the poorly-constrained density 
dependence of the symmetry energy |321 140| . Tuning this parameter has served to uncover correlations between the 
neutron skin of a heavy nucleus — such as ^"^Pb — and a host of both laboratory and astrophysical observables. 

A consistent mean-field plus RPA (MF-I-RPA) approach to the nuclear response starts with the calculation of 
various ground-state properties. This procedure is implemented by solving self-consistently the appropriate set of 
mean-field (i.e., Kohn-Sham) equations generated by the Lagrangian density given above [33J. For the various meson 
fields one must solve Klein-Gordon equations with the appropriate baryon densities appearing as their source terms. 
These baryon densities are computed from the nucleon orbitals that are, in turn, obtained from solving the one-body 
Dirac equation in the presence of scalar and time-like vector potentials. This procedure must then be repeated until 
self-consistency is achieved. What emerges from such a self-consistent procedure is a set of single-particle energies, a 
corresponding set of Dirac orbitals, scalar and time-like vector mean-field potentials, and ground-state densities. A 
detailed implementation of this procedure may be found in Ref. |41j . Having computed various ground-state properties 
one is now in a position to compute the linear response of the mean-field ground state to a variety of probes. In the 
present case we are interested in computing the isoscalar monopole response as probed, for example, in a-scattering 
experiments [37]. Although the MF-I-RPA calculations carried out here follow closely the formalism developed in 
much greater detail in Ref. [T3], some essential details are repeated here for completeness. 

The distribution of isoscalar monopole strength may be extracted from the imaginary part of a suitable polarization 
tensor that we compute in a consistent relativistic RPA approximation [42 l l43 j . That is, 

^L(9,c^) = ^|(vi/„|/5(q)|*o) ',5(l. -..„) = - -5 n™pA(q,q;^), (2) 

n 

where ^o is the nuclear ground state and ^„ is an excited state with excitation energy ujn — En — Eg. To excite 
isoscalar monopole modes a simple transition operator of the following form may be used: 

p(q)= /dV7/;(r)e-*i"-7V(r) ■ (3) 



Here p(q) is the Fourier transform of the isoscalar (baryon) density and 7" = diag(l, 1, — 1, — 1) is the zeroth (or 
timelike) component of the Dirac matrices. Such a transition operator is capable of exciting all natural-parity states, 
including the isoscalar monopole {EO) mode of interest to this work. In the particular case of an EO {J^ — 0+) 
excitation, the effective transition operator reduces to the following simple form: 

p,M)= /dV^(r)jo(gr)7V(r), (4) 



where Jq{x) — sin(a;)/a; is the spherical Bessel function of order zero. Finally, in the long-wavelength limit, the 
distribution of isoscalar monopole strength R{lu\ EO) may be directly extracted from the longitudinal response. That 

is, 

i?(c^;£;0) = V|(*„|r2|*o) ^S{uj-cj„) = lim f^] SL{q,uj;EO) . (5) 



Connecting the nuclear response to the polarization tensor is highly appealing as one can bring to bear the full power of 
the many-body formalism into the calculation of observables than can be directly extracted from experiment |42l I43j . 
Moreover, the spectral content of the polarization tensor is both illuminating and physically intuitive. However, 
enforcing the self-consistency of the formalism, while essential, is highly non-trivial. Yet maintaining self-consistency — 
in particular by using a residual particle-hole interaction that is identical to the one used to generate the mean- 
field ground state — is essential for the preservation of important symmetries and the decoupling of various spurious 
modes |13l 144) . Finally, given the critical role that certain moments of the distribution of strength play in our 
discussion, we close this section with a few essential definitions and relations. 

In general, the moments of the distribution of isoscalar-monopole strength are defined as follows: 
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In particular, the EWSR for the isoscalar monopole mode is given by [TS] 
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where M is the nucleon mass and p(r) is the ground-state baryon density normalized to the baryon number A. In 
essence, the power of the sum rule is that it relates a moment of the full RPA distribution to the mean-square value 
of the ground-state density. We note that the above classical sum rule is only valid in the absence of exchange 
and momentum-dependent forces [15] . Such forces modify the classical sum rules and their impact is traditionally 
accounted for by multiplying the right-hand side of Eq. (I?]) by a correction factor. Perhaps the best known case for the 
need of such a correction factor (k) is in the context of the photo-absorption cross section and the model-independent 
Thomas- Reiche-Kuhn sum rule (with ktrk ~ 0.2) |15| . As mentioned earlier, the classical sum rules were derived 
using a non-relativistic formalism so one may also need to correct for relativistic effects. We assume here that such 
relativistic effects may also be subsumed into such a correction factor. 

In addition to the EWSR a moment of critical importance to the present work is the inverse EWSR. This is because 
such a moment may also be computed from the ground-state density of the system. Indeed, 
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A=0 



However, the ground-state density in this case must be obtained by supplementing the mean-field potential with a 
"constrained" one-body term of the form Xr^. This can be easily implemented in the relativistic approach by adding 
such a harmonic one-body potential to the repulsive vector interaction. The addition of such a term shifts the weight 
of the single-particle orbitals to the interior, thereby leading to a more compact system. The inverse EWSR measures 
the (negative) slope of the mean-square radius at A = 0. Finally, the constrained energy is defined as: 
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Note that we have reserved the term "centroid energy" to the more conventional ratio of Eccn — Tni/mQ. 



III. RESULTS 



We start the section by displaying in Fig.[T] the distribution of isoscalar monopole strength for ^'^^Pb using the 
accurately calibrated FSUGold ( "FSU" ) parametrization [45] . Model parameters for this and the NL3 model are listed 
in Table |I] Note that the left-hand figure displays the energy weighted monopole strength whereas the right-hand 
panel the inverse energy weighted strength. We should also mention that due to the non-spectral character of the RPA 
approach [13] , the particle-escape width is computed exactly within the model. Given that the distribution of monopole 
strength is defined as the long- wavelength limit of the longitudinal response [see Eq. ([S])], we display the strength 
function for two small values of the momentum transfer to ensure the convergence of our results. Finally, the two insets 
display the cumulative sums relative to their corresponding sum rules computed from the constrained RMF approach as 
indicated in Eqs. (l7| and ((sl). The insets indicate that the RPA response accounts for about 90% of the corresponding 
sum rules. However, the constrained energies agree to better than 2%, namely, £'con(CRMF) = 13.50 MeV and 
-E'con(RPA) = 13.76 MeV. Note, however, that this 2% difference represents the largest discrepancy obtained in the 
present work (see Table [h]). 



2.0 



*Pb (FSU) 



q=0.095 fm"' . 
q=0.122fm-' - 




10 20 30 

(0 [MeV] 



*Pb (FSU) 




— Eeo„=13.76MeV 

— Eeo„=13.75MeV 



I I I 



10 20 30 40 i 
CO [MeV] 



^ ' ' ' 



40 



10 20 30 

CO [MeV] 



' ' ' ' 



1.0 



0.8 



> 






O 



-0.4 



^ 



-0.2 



40 



J3.0 



FIG. 1: (Color online) The energy weighted monopole strength (a) and the inverse energy weighted monopole strength (b) in 
■^°*Pb for two (small) values of the momentum transfer. The insets display the cumulative sums relative to the corresponding 
sum rules obtained from the constrained RMF approach [see Eqs. n\ and ([8|]. 
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TABLE I: Parameter sets for the two accurately calibrated relativistic mean-field models used in the text: NL3 |46l I47| and 
FSUGold [3S]. The parameter k, and the meson masses rris, m^, and nip are all given in MeV. The nucleon mass has been fixed 
at M = 939 MeV in both models. 



Before proceeding any further with the discussion of our results we should mention some subtleties associated with the 
extraction of the constrained energy from the RPA results. Essentially, the potential problems emerge from the lack of 
convergence of the integrals over the excitation energy lo defining both mi and m_i. Indeed, the lack of convergence 
of the energy-weighted and inverse energy-weighted sums can be clearly appreciated in the insets of Figs. fUa) and 
flTb), respectively. Note that this situation is not exclusive to the monopole resonance but extends to all excitation 
modes. To elucidate the problem we resort to an ideal Lorentzian distribution of strength that provides an accurate 
representation of the monopole strength in heavy nuclei, such as ^°®Pb. The Lorentzian distribution is defined as 
follows: 



R{uj) = Too 
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where toq is the unweighted-energy sum and F is the full width at half maximum. Note that the full unweighted-energy 
sum Too can only be recovered by extending the integral over the unphysical oj < region. That is, 



R{uj) du) = TO-o 
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However, by limiting the integral to the physical w > region, one may still account for most of the unweighted-energy 
sum — especially in the case that F^wq. Indeed, 
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(12) 



In the particular case of ^°^Pb with the parameters predicted by the FSUGold parametrization (a;o = 13.82 MeV and 
r = 2.29MeV) one can still account for about 97% of the unweighted-energy sum. However, the situation is radically 
different with the rrii and m_i moments as both integrals diverge; mi displays an ultraviolet divergence whereas r7i_i 
an infrared divergence — both logarithmic. Yet both integrals are well behaved if the integration region is allowed to 
be extended to the unphysical region. That is, 
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Note that in the case of an RPA distribution of an exact Lorentzian shape the centroid and constrained energies will 
be given by the following simple expressions: 
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(14) 



So how does one extract the constrained energy from a distribution of strength that is physically meaningful only 
for positive values of uj7 Given the previous discussion, we suggest to rely heavily on the unweighted distribution 
of strength as it is both well behaved and accounts for almost 100% of the sum rule (at least in the Lorentzian 
approximation). Thus, we select the upper limit of integration (wmax) in such a way that most of the integrated 
strength ttiq has been accounted for. For the case of the energy-weighted sum, the remaining contribution (from 
'^max to cx)) is assumcd to be exactly cancelled by the contribution in the unphysical region. Finally, to remove the 
w = singularity in to_i, we replace the low-energy tail by a rapidly falling exponential distribution in the interval 
0< a;< Wmin. Given that the divergence at both low and high excitation energy is logarithmic in nature, we found 
our results stable against small changes in both Wmjn and lUJmax- 
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TABLE II; Giant monopole resonance constrained and centroid energies for a variety of nuclei as predicted by the NL3 
parametrization [461 147] . Experimental data extracted from Refs. [371 1481 149] . All quantities are given in MeV. 
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TABLE III: Giant monopole resonance constrained and centroid energies for a variety of nuclei as predicted by the FSUGold 
parametrization [45; . Experimental data extracted from Refs. J37l l48l I49j . All quantities are given in MeV. 



Having explained how we extract RPA centroid and constrained energies, we now continue with the discussion of our 
results. Although the use of mean-field methods for light nuclei may be questionable — especially since the distribution 
of monopole strength may be strongly fragmented — we have performed self-consistent RPA calculations for nuclei 
ranging from ^^O to ^^^Pb using the NL3 and FSUGold energy density functionals. Gonsidering these two functionals 
is useful because, although they both provide accurate descriptions of a variety of nuclear ground-state properties, their 
predictions for infinite nuclear matter are significantly different pi]. In particular, the incompressibility coefficient of 
symmetric nuclear matter predicted by NL3 is about 271 MeV whereas that of FSUGold is 230 MeV. We provide our 
results in both tabular form in Tables[lT| and |IIH and graphical form in Fig. [2] First, we note that the NL3 results for 
the centroid energies are systematically higher than the corresponding FSUGold results, as expected given the higher 
incompressibility coefficient of the NL3 interaction. The NL3 centroid energies also are systematically higher than 



the experimental data except in the case of ^"^^Pb. The case of ^°^Pb is unique because the large incompressibility of 
NL3 is compensated by a similarly large slope of the symmetry energy; this accounts for the accurate prediction of 
the centroid energy in ^°*Pb — but not in ^°Zr [TU]. In essence, the FSUGold interaction was conceived with the goal 
of better constraining the incompressibility of neutron-rich matter by incorporating information on giant resonances 
into the calibration procedure. However, in the present work we aim to go one step further by exploring the feasibility 
of incorporating monopole energies directly into the calibration of the quality function. The results presented in both 
of the tables are extremely encouraging as they suggest differences of at most 2% between energies computed with 
the constrained RMF approach and those extracted from the RPA strength function. For reasons that at present we 
do not understand, the discrepancy between constrained and RPA results is systematically higher in the case of the 
FSUGold parametrization. Nevertheless, the accuracy of our results gives us confidence that in the future, constrained 
GMR energies for nuclei ranging from ^°Zr to ^°^Pb — and including a variety of Sn-isotopes [T^J [SOI [SI] — may be 
directly included in the accurate calibration of energy density functionals. Such a procedure will not only be able 
to better constrain the density dependence of the equation of state around saturation density, but may be helpful in 
partly removing the reliance on "pseudo-data" (i.e., on bulk properties of infinite nuclear matter). 
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FIG. 2; (color online) Centroid and constrained energies for a variety of nuclei as predicted by the NL3 
models. Experimental centroid energies are taken from Refs. [371 1481 149] . 
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IV. CONCLUSIONS 



The sum-rule approach provides a powerful tool for the analysis of nuclear excitation spectra [15] . Particularly useful 
are the energy-weighted sum rules — with the TRK sum rule as its best known exponent — as they can be determined 
largely model-independently from a few well-known ground-state properties. Moreover, additional sum rules, such as 
the inverse energy weighted sum rule employed in this work, may be combined to estimate the mean excitation energy 
of the resonances from intrinsic ground-state properties. However, whereas the inverse EWSR may be extracted from 
a constrained calculation that relies exclusively on first-order perturbation theory, the model-independent EWSR is 
based on commutation relations between an appropriate transition operator and a non-relativistic Hamiltonian. Thus, 
the sum-rule approach has been limited to the non-relativistic domain; for a recent and successful implementation of 



such techniques see Ref. fTO and references therein. Indeed, in Ref. TO| highly accurate Skyrme-based mean-field plus 
RPA calculations are successfully compared against constrained predictions for the mean excitation energy of various 
resonances. 

In this contribution we suggested that the constrained approach should also be valid in the relativistic domain. We 
justified such a conjecture by invoking the basic tenets behind density functional theory. DFT asserts that the exact 
ground-state energy and density may be obtained through the minimization of a suitable energy density functional. 
Moreover, within DFT, the linear response of the system may be computed from the small density fluctuations around 
the variational minimum. In particular, no demands are made on whether the functional should be of a relativistic or 
non-relativistic character. Indeed, as long as the energy density functional is accurately calibrated to physical data, 
the precise form of the functional becomes a matter of convenience. We believe that assessing the degree of agreement 
in the relativistic domain between the predictions of the constrained approach and the RPA results is reason enough 
to implement the present study. However, our study also has an important practical component. Given that we aim 
to eventually supplement the calibration of energy density functionals with experimental data on nuclear collective 
modes, the validity of the constrained approach becomes essential, as it is impractical to incorporate this information 
from a detailed RPA calculation. 

Relativistic mean-field plus RPA calculations were performed for the distribution of isoscalar monopole strength 
for magic or semi-magic nuclei ranging from ^^O to ^''^Pb using two accurately calibrated density functionals (NL3 
and FSUGold) . Predictions were made for both centroid and constrained energies from the various moments of the 
distribution. RPA calculations of this kind are computationally expensive as the strength distribution is finely mapped 
over a wide energy range. In contrast, the extraction of monopole energies from the constrained approach is numerically 
expedient as it relies exclusively on the self-consistent evaluation of ground-state properties. As previously conjectured, 
we find excellent agreement between the two approaches. Indeed, for all nuclei under consideration the discrepancy 
between the two approaches amounts to less than 2%, a quantity that is comparable to current experimental errors. 

In summary, motivated by ideas from density functional theory we proposed an extension of the constrained ap- 
proach to the relativistic domain. Relativistic constrained calculations of monopole energies were favorably compared 
against the corresponding predictions from the relativistic RPA approach. These results are highly encouraging and 
show the promise of incorporating nuclear excitations in the calibration of future energy density functionals without 
incurring in an unaffordable computational cost. 
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